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ABSTRACT 

We present limits on the molecular gas content of Lyman Break Galaxies (LBGs) 
at z ~ 5 from observations targetting redshifted CO(l-O) and CO(2-l) line emission. 
We observed a single field containing eight spectroscopically-confirmed z ~ 5 LBGs, 
seven of which are contained within a narrow (z=4.95 ± 0.08) redshift range and the 
eighth is at z — 5.2. No source was individually detected. Assuming the CO to H2 
conversion factor for vigorous starbursts, we place upper limits on the molecular gas 
content of individual z ~ 5 LBGs of M(Ha) < 1O 1O M0. From a stacking analysis 
combining all of the non-detections, the typical z ~ 5 LBG has an H2 mass limit com- 
parable to their stellar mass, < 3.1 x 10 9 M Q . This limit implies that, given the star 
formation rates of these systems (measured from their UV emission) , star formation 
could be sustained for at most ~ 100 Myr, similar to the typical ages of their stellar 
populations. The lack of a substantially larger reservoir of cold gas argues against the 
LBGs being UV luminous super starbursts embedded in much larger UV-dark systems 
and as a result increases the likelihood that at least those LBGs with multiple compo- 
nents are starbursts triggered by mergers. The sources responsible for reionization are 
expected to be starbursts similar to these systems, but with lower luminosities, masses 
and consequently with star formation timescales far shorter than the recombination 
timescale. If so, the ionized bubbles expected in the IGM during the reionization era 
will only infrequently have UV-luminous sources at their centres. 

Key words: galaxies: high-redshift - galaxies: starburst - galaxies: star formation - 
radio lines: galaxies 



1 INTRODUCTION 

Detailed observations of the earliest galaxies are necessary 
if we are to form a complete picture of galaxy formation and 
evolution. While increasing numbers of spectroscopically- 
confirmed galaxies are being discovered at z ~ 5 and above 
(within ~lGyr of the Big Bang), almost all are discovered 
through rest-frame UV emission originating from strong on- 
going star formation. Unfortunately, this tells us little about 
the (potentially dominant) UV-dark baryonic component of 
these galaxies and consequently limits our understanding of 
star formation in the high-z universe. 

Lyman-break Galaxies (LBGs) form a s ubstantial part 
of cu r rent samples of z > 5 ga laxies (e.g. IVanzella et al.l 
120091 : iDouglas et all 120091 . |201Ch . They are identified via 
their UV continuum emission longward of 1216A in the 
rest-frame, which arises from hot young stars formed in un- 
obscured starburst regions. While Spitzer-based follow-up 
observations have made some progress in exploring and plac- 
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ing limits on their older underlying s tellar populations (e.g. 
lEvles et al]|2007l : IVerma et al-lhoOTT ). we know little about 
their gas content, which is a crucial diagnostic of the du- 
ration of the ongoing starburst and of the nature of the 
LBGs themselves. The galaxies have unobscured star forma- 
tion rates of a few x lOMQyr -1 arising from regions with 
a typical surface area of ~ lkpc 2 . The strong wind gen- 
erated by such a starburst can potentially limit the avail- 
able fuel for continuing star formation. Clearly, the amount 
of gas available to the ongoing starburst is key to an un- 
derstanding of the nature of the burst and the future of 
star formation act ivity within the system. HST imaging of 
z ~ 5 LBGs (e.g. IDouglas et al.l [2OI0I ; IConselice fc Arnold! 
120091 ) shows that many sytems have multiple UV compo- 
nents and extended, distorted morphologies on scales of a 
kpc or more. With the available optical and near-IR data, it 
is currently impossible to determine whether such structures 
imply that the LBGs originate in mergers, or are individual 
super starburst regions embedded in much larger UV-dark 
systems, as found in low r edshift Lyman break analogues 
(Overzier et all l|2008l ). see IDouglas et al.1 l|201Ch for a dis- 
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cussion). If the latter scenario is correct, then one would 
expect considerably more cool and cold gas present in the 
immediate environment of the LBGs than is present in the 
stellar mass produced by the ongoing starbursts. 

In this paper we probe the cool gas component of a sam- 
ple of z ~ 5 L BGs drawn from the ESQ Remot e Galaxy Sur- 
vey (ERGSEI, iDoudas et ai]|2007t 120091 I2010T I. Two of the 
ten 40 arcmin 2 ERGS pointings display a large over-density 
of spectroscopically-confirmed UV bright sources over nar- 
row (Az ~ 0.1) redshift ran ges. The first of thes e was the 
subject of a previous letter IjStanwav et alj 1200a ) in which 
we discussed the identification of a UV-dark molecular line 
emitter at z ~ 5. Here, we target one field, J1054. 4-1245, 
which contains many spectroscopically-confirmed LBGs in 
a two-arcminute diameter region. As in our previous work, 
we target the field using the Australia Telescope Compact 
Array (ATCA), but this time use the new Compact Array 
Broadband Backend (CABB) which probes a Az ~ 0.6 red- 
shift range in a single exposure (at 38GHz). This is approx- 
imately 15 times broader than in our first work and easily 
encompasses the entire velocity range probed by the LBG 
over-density. 

Throughout this paper, unless otherwise stated, all 
magnitudes are on the AB scale, and the cosmology used 
is Fo=70kms _1 Mpc" 1 , Qa = 0.7 and M =0.3. 
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Figure 1. The spatial position of spectroscopically-confirmed 
LBGs in our target field. Red diamonds indicate LBGs in the 
narrowest redshift over-density, while the black cross shows the 
eighth LBG in this region at slightly higher redshift. Our ATCA 
pointings are indicated by the 7mm (solid circles for 2009 obser- 
vations and dot-dashed blue circles for 2010 observations) and 
12mm (dashed circle) half power beam widths. These encompass 
all eight LBGs. 



2 OBSERVATIONS 

Observations were carried out with the ATCA over two sep- 
arate runs, one in May 2009 as part of project C1954 and 
another in March 2010 as part of project C2297. The pri- 
mary goal of the project was to search for UV-dark molec- 
ular line emitters across our target field and the results of 
that study will be presented in a future paper. The array 
was in the compact H168 configuration, using both North- 
South and East- West baselines. We utilized the new CABB 
correlator with two intermediate frequency (IF) bands, each 
with 2GHz/2048 channel configurations. We targeted CO(2- 
1) and CO(1-0) transitions at the LBG over-density redshift. 
For CO(2-l) transitions we tuned the IF bands to 36.70 
GHz and 38.72 GHz to allow a simultaneous survey range of 
4.81 < z < 5.44. We bin six adjacent 1 MHz channels to in- 
crease signal to noise, and we obtain a spectral resolution of 
~47km s _1 . During the first run, observations were taken in 
six 8 hour periods between 2009 May 1 and 7. Three point- 
ings were observed in order to target the maximum possible 
number of LBGs in the field. During the second run two ad- 
ditional pointings were observed in four 8 hour periods be- 
tween 2010 March 21 and 23 (see figure[l|). A nearby bright 
source, PKS1054-188, was observed every 15 minutes to de- 
termine the phase stability and the pointing accuracy was 
was checked every hour. Primary flux calibration was carried 
out on the standard ATCA calibration source, PKS1934- 
638, each night. The half-power-beam-width (HPBW) of the 
ATCA at ~37 GHz is 74" and the restoring beam, for nat- 
ural weighting and this configuration, is 7.3" x 4.8". For 
CO(1-0) transitions we tuned the IF bands to 19.12 GHz and 
21.12 GHz giving redshift coverage of 4.23 < z < 5.34. We 
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observed a single pointing for two nights on 2009 May 8 to 9, 
encompassing all five ~37GHz pointings (see Figure [TJ. As 
before, we bin three adjacent 1 MHz channels to increase sig- 
nal to noise which gives a resolution of ~47km s" 1 (match- 
ing the velocity resolution of our high frequency data), use 
nearby source 1054-188 for secondary flux calibration and 
PKS1934-638 for primary flux calibration. The HPBW at 
~20 GHz is 2.5' and the restoring beam is 15.9" x 10.9". At 
z ~5 this corresponds to a beam size of > 50kpc. 

We expect any detectable line emission from our sources 
to have velocity width ~150km s _1 (the only two known 
high redshift, non AGN, CO line- emitting galax i es at z ~ 5 



have line widt hs of ~ 160km s 1 , Coppin et al.l (J2010J) and 



~ 110km s _1 , IStanwav et al.l 120081 ) and we require detec- 
tions to be statistically significant (> 2<r) in each of 3 adja- 
cent channels. Therefore by binning our data to 3MHz and 
and 6MHz at the lower and higher frequency settings, a de- 
tection in three adjacent bins will correspond to a ~150km 
s _1 -wide line. 

Total integration times were 16 hours at each frequency 
giving an rms noise of ~0.11 mjy/beam at 19 GHz and 
~0.17mJy/beam at 37GHz in each 47km s _1 channel. 



3 LIMITS ON LBG CO EMISSION 

Our observations target eight spectroscopically-confirmed 
z ~ 5 LBGs which are spatially distributed within the half 
power beam widths of individual pointings at both ~37GHz 
and ~19GHz. Seven lie in the 3-dimensional over-density 
(z ~ 4.95 ± 0.08) of LBGs in the larger ERGS field and 
the eighth is at slightly higher redshift (z ~ 5.2). None of 
these galaxies show detectable line or continuum emission 
at either of the redshifted CO transitions targeted. For the 
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Velocity offset from Ly-alpha redshiit, km/s 

Figure 2. ID spectrum of a typical LBG in our target region centered on the rcdshifted CO(2-l) line with redshifts assuming the Ly-a 
redshift. The blue solid line is the target spectrum (binned to 6MHz channels) with errors displaying the la rms error in a 20" X 20" 
region around the LBG position. The red dashed line shows the target spectrum convolved with a 150km s" 1 (FWHM) Gaussian. 
The shading highlights the region in which CO lines will found given typical offsets between Ly-a and systematic redshifts. The most 
significant positive peak in the above spectrum at +650km s _1 is still no more than 2a away from zero and consequently, given the 
number of independent bins displayed, is entirely consistent with noise. 



line emission, each extracted 1-dimensional spectrum was 
examined for emission at the frequency expected for red- 
shifted CO emission (using the redshift determined from ei- 
ther Lya em ission or the break a t Lya detected in optical 
spectroscopy iDouglas et al.l 120101). However, in a study of 
z ~ 3 star forming galaxies, ISteidel et al.l (|20ld ) show that 
the Lya redshift can differ from the systematic redshift de- 
rived from interstellar absorption lines and other emission 
lines by up to ~ 600kms _1 , with the Lya redshift being 
systematically higher. Consequently, we additionally deter- 
mined limits to any line flux at frequencies corresponding 
to velocity offsets of up to 600kms _1 shortward of the Lya- 
derived value (see figure [2]). 

To search for CO line emission, 47 kms - -wide slices 
of the data cube were extracted and the root-mean-square 
variation was ascertained in each for a region of ~20 x 20 
arcsec box centred on the position of each LBG. Any real 
emission line is likely to be broader than 50 km s _1 as noted 
above. Consequently, we define a detection as three consec- 
utive 47 kms - spectral bins at least 2a above zero flux 
in the region between and 600 kms -1 shortward of the 
Lya redshift. Using these criteria, no individual source was 
detected. To determine a characteristic line-limit for each 
source, eleven values of the rms over three consecutive bins 
in the same to 600 kms -1 region were obtained. As the 
variation between these was small, we took the average 
value of these as characteristic of a given source. This re- 
sults in a typical flux limit of < 60 mjy km s _1 for the 
CO(2-l) lines and < 40 mjy km s -1 for the CO(1-0) line. 
These correspond to luminosity limits in each line of typ- 



ically 1.1 x 10 1U K km s 



1 PC 2 



at z ~ 4. 95 for CO(2-l) 



and 2.6 x 10 10 K km s" 1 pc 2 for CO(1-0) (|Solomon et all 
1 19971 ). Converting this to a limit on the molecular gas mass 
requires an appropriate conversion factor. In the absence 
of a directly measured value at the highest redshifts, we 
use the commonly-used conversion factor derived from lo- 



cal st rongly star forming galaxies (|Solomon fc Vanden Bouti 
120051 ) of M gas /L' co = 0.8M Q (Kkm s -1 pc 2 ) -1 given the 
strength of the starbursts in these LBGs. Assuming that 
the emitting medium is optically thick and thermalized, 
the line lumin osity is independent of transition [i.e. L co (2- 
1)=L' CO (1-0), ISol'omon fc Vanden Boutl (|2005l )). The typi- 
cal CO(2-l) limit places a constraint on the gas mass of 
Mh 2 % 8.9 x 10 9 Mq . This is comparabl e to the stellar mass 
content of LBGs at this redshift {e.g. IVerma et al.l 120071 ; 
IStark et alll2009h . The constraint placed by the CO (1-0) 
line is a factor of two to four times higher depending upon 
the source. 

Of course, if the line width is much narrower than we 
assumed, we may have missed emission using these criteria. 
However, we note that the above flux and luminosity lim- 
its are appropriate for a 6a non-detection of any emission 
confined to a single channel. In fact, no channel in the 600 
kms -1 range in any of the spectra deviates from zero by 
more than 2.5a, so the limit quoted above is robust. 

In addition to determining limits on individual sources, 
we can derive deeper limits on the "typical" z ~ 5 LBG 
by creating an average stack from the spectra of the eight 
individual sources, having shifted each spectrum to the same 
effective redshift. Concentrating on the CO(2-l) data and 
using the same criteria as for an individual spectrum, no 
line was detected in the averaged spectrum. The limits on 
the flux, line luminosity and H2 mass of a "typical" LBG 
are tightened to 5*coA^ < 17.5 mjy km s , £00(2-1) < 
3.9 x 10 9 K km s -1 pc 2 and M H , < 3.1 x 10 9 M at z = 4.95 
(2a). 

However, this process is not completely straightforward 
as the true redshifts for the sources may not exactly match 
the Lya-derived values, as noted above. The offsets between 
the Lya and true values may vary by ~ 300km s -1 between 
objects. Consequently, any combination of weak lines may 
be averaged in such a way that they are smeared out in 
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Table 1. Properties of LBGs in target field (Positions given in J2000). NOTES - 1 Average of the total of the la rms errors in 3 adjacent 
bins (total in lSOkms^ 1 width) for all bins, 600kms _1 blue-ward of the CO(2-l) line position if optical redshifts are correct. 2 Limit of 
CO luminosities derived from rms errors for an unresolved ~150km s —1 line (three channel s at la rms limit). ^Inferred gas mass derived 
from conversion factor for local infrared-luminous galaxies (Sol omon &: Vanden Boutll2005h . 
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Figure 3. Histogram of the distribution of flux limits from 823543 
possible iterations of stacking LBG spectra. Eight LBG spec- 
tra were combined varying the relative offset in redshift between 
seven of the spectra in six steps, +150, +100, +50, 0, -50, -100 
and -150 kms -1 relative to the eighth spectrum. The dashed line 
shows the flux limit for the eight spectra stacked at the optically- 
derived redshifts. The distribution of flux limits is Gaussian- 
distributed as expected given the noise characteristics of indi- 
vidual velocity channels. 



frequency and do not reinforce each other to become a de- 
tection in the average spectrum. To explore this we made 
multiple average spectra, varying the relative offset in red- 
shift between each spectrum in seven steps, +150, +100, 
+50, 0, -50, -100 and -150 kms -1 , in total 7 7 or 823543 av- 
erage spectra. We then take the ten flux limits in the -600 to 
0km s _1 region from each spectrum. The distribution of flux 
limits for these spectra are shown in figure[3] As can be seen, 
using the Lyev redshifts does not significantly underestimate 
the true value of the limit. 

For completeness, we carried out the same procedure, 
but looking at features contained within a single bin as 
would be the case for extremely narrow lines. The highest 
flux in one of the -600 to kms - bins in the average spectra 
was 12.6 mjy kms -1 (consistent with the noise character- 
istics given the 823543 different realisations), considerably 



less than the above limit of S'coAi' < 17.5 mjy kms -1 and 
so we consider that limit to be robust even for a very narrow 
line. 



4 DISCUSSION 

Assuming that the conv ersion between CO lumin o sity a nd 
H2 mass determined by ISolomon fc Vanden Bouj i|2005l ) is 
applicable (the intensity of the starbursts in these LBG in- 
dicates that it is), the non-detection of CO emission from 
these LBGs places interesting constraints on their nature. 
The similarity between the typical stellar masses of such 
system s as determined fro m multiband photometry of 2 ~ 5 
LBGs (|Verma et al.ll2007l ). and the limit to their typical H2 
mass (a few times larger) would indicate that the systems are 
observed about half way through the ongoing burst of star 
formation assuming that the current star formation rate is 
maintained and the conversions from gas to stars approaches 
~ 100 per cent efficiency. However, the efficiency with which 
gas is converted into stars is typically no mor e than ten per 
cent, and potentially only one or two per cent (JLehnert et al.l 
2009), thus suggesting that these sources are much more 
than halfway through their life cycle. This is unsurprising 
as timescales of ~ 10 Myrs will be required to produce and 
sustain the observed UV continuum fluxes from a popula- 
tion of O and B stars. Dividing the limiting gas mas s by the 
typical star formation rate from lVerma et al.l (|20070 gives a 
timescale of ~ lOOMyr , with the typ i cal ag e of a stellar pop- 
ulation determined bv lVerma et al.l (120071 ) being a few tens 
of Myr. As this implicitly assumes an unrealistic complete 
transformation of gas into stars, even if the CO luminos- 
ity to H2 conversion factor is higher than the Solomon & 
Vanden Bout value, this is a robust upper limit to the star 
formation timescale at this level of activity. 

The apparent lack of a comparatively large reservoir 
of molecular gas in and around the UV-luminous system 
strongly argues against the LBGs being unobscured super 
starbursts embedded in much more massive and extended 
underlying systems. Any such system should contain suffi- 
ciently large amounts of enriched molecular gas to be de- 
tectable here. Consequently, it provides support for the al- 
ternative hypothesis of a merger origin for at least those sys- 
tems with multiple spatially-distributed UV luminous com- 
ponents. If all z ~ 5 LBGs are triggered by mergers with 
other galaxies, our lack of a single detection implies that few 
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involve comparatively massive systems with substantial ob- 
scured star formation or large amounts of processed molec- 
ular gas. Additionally, the most active phases of a merger 
occur when the gas has been consumed and is concentrated 
in the central regions of the systems. If LBGs are triggered 
by mergers, and given the need to build UV continuum fluxes 
through star formation, this may take substantial time. 
Therefore in this model LBGs may only be detectable (in 
the rest frame UV) towards the end of their starbursts. Ob- 
served luminous starburst galaxies at high- 2 typically have 
short gas depletion times relative to their ages, similar to 
that seen in these sources. The lack of detectable 850 micron 
continuum flux in APEX/LABO CA observations of differ- 
ent but essentially similar LBGs (jStanwav et al. 20101) gives 
furthe r support to this model. As noted in IStanwav et al.l 
(|2010r) the conversions from continuum flux limits to star 
formation rates are currently ambiguous. However, the lack 
of detectable dust emission from these similar sources lends 
support to the idea that they do not have large reservoirs 
of cool and cold material, thus strengthening the idea that 
these sources are not embedded in more massive obscured 
systems. 

There has been re cent discussion 

llPapadopoulos fc Pelupessvl |2010| ; iGnedin fc Kravtsovl 
[2010J) of whether there should be deviation away from the 
low-redshift Schmidt- Kennicutt (SK) relation (ISchmidtl 
1 1953 ; iKennicut j 1 1998T ) at high redshift. Here we note that 
for unobscured star formation rates of a few tens of solar 
masses per yea r arising from reg i ons of approximately 
~ lkpc (e.g IVermaetal.1 120071 : iBouwens et~aI1 12004 
iBremer et al.l 120041 ). a limit of Mh 2 < 3.1 x 10 9 M in the 



same area is consistent with the low-redshift SK relation. 
The synthesised ATCA beam at 38GHz probes an area of 
well over 100 kpc 2 centred on each LBG and so if there 
was an underlying larger UV-dark galaxy within which 
each LBG was embedded, any significant obscured star 
formation in that galaxy would lead to deviation from the 
nearby SK relation. As simulations imply that the relation 
gets steeper with redshift (more gas per unit star formation, 
e.g. IGnedin fc Kravtsovl [2010J), the deviation from the 
SK relation would only increase if the nearby relation is 
inappropriate. 

Given the comparatively brief elapsed time between the 
end of reionization and the redshift explored here, the lim- 
ited gas reservoirs available for star formation in these sys- 
tems and their relatively short star formation timescales 
have an important consequence for the reionization pro- 
cess. Reionization is likely to be dominated by as yet 
undetected faint, low mass UV-luminous star bursts (e.g. 
iLehnert fc Bremerl 120031 ; ISalvaterra et al.ll2010l '). As a con- 
sequence of their lower mass, these sources will have even 
shorter star formation lifetimes than the LBGs observed 
here. At all redshifts the recombination time for the IGM 
is longer than the Hubble time and hence also much longer 
than the UV-luminous lifetime of a low mass starburst. Con- 
sequently, most of the ionized bubbles expected to exist in 
the otherwise neutral IGM during reionization will only oc- 
casionally have UV-luminous sources at their centres, and 
most starbursts giving rise to the ionising flux will have long- 
since ceased to actively form stars and therefore have faded 
in the UV. 
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